YCTP-N1-95 
January, 1995 



Quantum Dissipation 



Aurel BULGAC a,fe '[|, Giu DO DANG fc f] and Dimitri KUSNEZOV c f| 

a Department of Physics, FM-15, University of Washington, Seattle, WA 98195, USAf\ 

b Laboratoire de Physique Theorique et Hautes Energies, 
Universite de Paris-Sud, Bat. 211, 91405 Orsay FRANCE 

c Center for Theoretical Physics, Sloane Physics Laboratory, 
Yale University, New Haven, CT 06511-8167, US A 



Abstract 

We address the question of the microscopic origin of dissipation in collective 
motion of a quantum many-body system in the framework of a parametric random 
matrix approach to the intrinsic dynamics. We show that the fluctuation-dissipation 
theorem is generally violated and, moreover, energy diffusion has a markedly non- 
Gaussian character and the corresponding distribution has very long tails. Such 
features do not support a Langevin or Fokker-Planck approach to dissipation in 
collective nuclear motion. 
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While theory and experiment have gone a long way in the study of the collective 
nature of large amplitude nuclear motion, the theoretical understanding of the coupling 
between the collective and the intrinsic degrees of freedom is still in its infancy. Most of 
the approaches are more or less phenomenological in nature. Hill and Wheeler JTJ sug- 
gested almost forty years ago that Landau-Zener transitions are at the origin of nuclear 
dissipation. Over the years there have been a relatively large number of studies of this 
particular mechanism |J and the range of results is equally diverse. Implicit in this inter- 
pretation is the presumption that irreversibility is of quantum origin. Even though there 
exist quantum approaches (the linear response model ||, the hopping model [|J, path in- 
tegral method |J and others), many formulations are basically classical , e.g. the so called 
"wall formula" 0. These include the more pragmatic phenomenological models, such as 
the Langevin equation or Fokker-Planck equation J7j, Maxwell's model for friction with 
memory effects [§, and to a certain extent kinetic approaches, e.g. two-body dissipation 
mechanism ||. In an analysis of a generic problem of coupled slow and fast degrees of 
freedom, Berry and Robbins f|] obtained friction for the slow subsystem only by treating 
the entire system classically, attributing it to "a clash between the essence of quantization, 
namely the discrete spectrum of frequencies, and the essence of chaos, namely mixing and 
a continuous spectrum extending to zero frequency" . The present status of our under- 
standing of the microscopic origin of dissipation in many-body quantum systems is thus 
rather unsatisfactory. In this letter we explore the nature of dissipation in a many-body 
system, using parametric random matrix theory, which allows for direct solution in many 
cases of the quantum dynamics of the system. 

The physical systems we explore are many-body systems which exhibit excitations 
on two distinct time scales, described by collective (slow) and intrinsic (fast) degrees of 
freedom. In order to address how energy is transferred from the slow (X, P) to the fast 
(x,p) modes, we will assume that the slow modes evolve classically at constant velocity 
Vq according to X(t) = Vot. (Although we do not consider it here, this restriction can 
be lifted, and the more general problem solved using the results and methods presented 
here.) As a consequence of this assumption, we can solve for the quantum dynamics of 
the fast subsystem, and even obtain analytic results for situations which are analogous 
to the conventional adiabatic and diabatic limits. The intrinsic system is defined by its 
matrix elements, and is taken as complex, described by its average level density, p(E), 
and its spectral fluctuations, in the form 

H(X) = H + H 1 (X). (1) 

Here H is chosen to be a diagonal N x N matrix, defining the average density of states, 
with (k\Ho\l) = [Hq]m = EkSki- In the basis of the eigenstates of H , we define Hi(X) as a 
parameter dependent, N x N real Gaussian random matrix, which is completely specified 
by its first two moments 

[H^XMH^Y)]^ = [SikSfl + SuSjdfijiX-Y), 

[Hi(X)] kl = 0. (2) 
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Tij{X — Y) is a "bell-shaped" correlation function with a characteristic width Xq, and the 
overline stands for the ensemble average. The dependence on i,j allows for the description 
of banded matrices, where an effective number of states N < N can be coupled by Hi(X). 
Such a parametrization implies that correlations between different instantaneous spectra 
corresponding to different 'shapes' X are effective only within a distance rs Xq. The 
average level density for each fixed shape X is given by Hq, while its spectral fluctuation 
properties (in this case GOE, but all our formalism applies to GUE as well) are determined 
by Hi(X). We use a convenient parametrization of the correlator JF^(X) introduced in 
Refs. [g 

X 



FiA X ) = i , . . . ex P 
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Here F(x) = F(-x) = F*(x) < 1, F(0) = 1 and W , k (iV ~ «op(e)) and X are 
characteristic to the given system. The instantaneous spectra of a Hamiltonian H(X) 
with constant average level density, [H Q ] k i = kdki, is shown in Fig. 1 for gaussian F(x) = 
exp(— x 2 /2) (top) and exponential F(x) = exp(— \x\) (bottom) correlations. Notice that 
a conventional adiabatic limit does not exist for the exponential, as the individual energy 
levels undergo Brownian motion, and are not smooth on any time scale. 

The time evolution of this system is found by solving the time-dependent Schrodinger 
equation: 



i[)(t) = Texp 



- % - jf dsH(X(s)) 



#)=WW(0). (4) 



where T is the time-ordering operator, and U(f) the propagator. (We assume that the 
initial state ^>(0) is uncorrelated with the Hamiltonian H(X(t)) at later times.) The 
average propagator U (t) —U(t), found by using Eqs. (2) and resumming all leading order 
diagrams in perturbation expansion of U(t) in the limit JVo > 1, satisfies the system of 
integral equations [11]: 

U k {t) = U 0h (t)--2 ds x ds 2 U k (s 2 )U k (t - st) 
n Jo Jo 

N 

x Y, Fkn(X( Sl ) - X(s 2 ))U n ( Sl - s 2 ), (5) 

n=l 

where U ^(r) = exp(—iekt/h). In order to compute averages of observables, we introduce 
the set of generalized occupation number probabilities 



fV k (t u t 2 ) = (ip(ti)\k)(k\ip(t 2 )) 



= Y(^(ti)\k)(k\U(t 2 )\l) ni (0), (6) 
i 

where riiit) = jVi(t,t) is the occupation probability of the state \l). Afk(ti,t 2 ) satisfy the 
following set of integral equations 

fVk(h,t 2 ) = U;(tx)U k (t 2 )n k (0) + -2 [ 1 ds 1 [ 2 ds 2 Y,M(si,S2)Fik(si- s 2 ) 

h Jo Jo i 

xUZ(t 1 -8 1 )U k (t 2 -s 2 ). (7) 
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These equations specify the time evolution of the system, and we will consider (z) the 
numerical solutions of (5)-(7) and the velocity dependence of the diffusion constant, and 
(ii) the extension of the formalism to the regime 1 <C N < N = oo, where we find 
analytic limits and a great simplification of the formalism as well. 

The first situation we study is that of constant average level density ([ifo]jw = k5ki), 
as in a stadium billiard. Eqs. (5)-(7) have been solved numerically for iV = 101 levels, a 
bandwidth N = 21, a gaussian correlation F(x) = exp(— x 2 /2) with X = 1, and initial 
conditions n^(0) = 6% 5%. The resulting occupation numbers n^(t) are shown in Fig. 2 
for the cases of fast (Vq = 4, top) and slow (Vq = 1/16, bottom) driving velocities. As 
the results are symmetric with respect to the index k, rikit) = n W 2-k(t), only k — 1 — 51 
are shown, counting from the bottom of the figure. One might expect that even a small 
driving velocity would result in a complicated time evolution, as the Hamiltonian is time- 
dependent, and has many small gaps in the instantaneous spectrum, where Landau- 
Zener transitions might occur and thus induce "irreversibility" 0. Actually, as we have 
discussed at length in Ref. jllj, this mechanism, which has been advocated in many 
previous treatments ||, is valid only for isolated level crossings and thus is unrealistic 
when there are many non-isolated ones as shown in Fig.l. 

In Fig. 2, one can clearly distinguish two time scales: a relatively rapid initial transient 
evolution, followed by a much slower one. While the initial transient behaviour is almost 
identical in both cases, governed by the same spreading width I>, the long time behaviour 
is strikingly different. For small driving velocities, the time evolution rapidly equilibrates, 
and can be understood in terms of the Vq — > limit, corresponding to constant random 
matrix theory |TT| . For large velocities there is a steady evolution to a different probability 
distribution. The initial transient behaviour arises only because our initial occupation 
probabilities rik(0) did not originate from an instantaneous eigenstate of H(0) (detailed 
discussion on initial conditions will be presented elsewhere ||10||). The subsequent long time 
behaviour is due to the explicit time dependence of the Hamiltonian H(t) and would be 
absent for a time independent one. 

The diffusion process associated with these time evolutions can be characterized by 
the energy variance, A E (i), and the energy diffusion constant D(Vq), defined by 



Ag(t) = (ip(t)\[H(t) — E(t)] 2 \ip(t)) « const +2D(V 2 )t. (8) 

t — ► 00 

In Fig. 3, D(Vq) can be seen to exhibit quadratic velocity dependence, in contradiction to 
previous claims 0. In the case we consider here of a symmetrical initial distribution rik(0) 
and constant level density, the average energy E(t) = (ip(t)\H(t)\tp(t)) is constant, hence 
the reaction force on the slow system, in particular the friction force, exactly vanishes. 
This is consistent with a fluctuation-dissipation theorem in the following sense. Expressed 
as 7 = (3D, where 7 is the "friction" coefficient and f3 = 1/T = d\ap{e)/de = is the 
inverse thermodynamic temperature, we have the expected result dE(t)/dt = 7 = 0. 
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We shall outline briefly a further extension and simplification of the formalism in the 
limit N — > oo, i.e. in the spirit of the standard constant random matrix theory, by 



introducing the characteristic functional [IT 

N(h,t 2 ,r) 



(V>(*i)|exp 
1 



tH (r-t 1 + t 2 ) 



h 



2-Khp(e k ) 
x exp 



J dTfS(t 1 ,t 2 ,T 

ie k {r -ti + t 2 ) 



Af(t,t, 



T 



cxp 



IT) 



Y,m)\m\m))h 



h n n\ 



(9) 
(10) 



where ((ip(t)\Ho \4>(t))) are cumulants. Af(ti,t 2 ,r) satisfies the evolution equation 



Af(t u t 2 ,T) = a*{t 1 )a{t 2 ) + 



1 


i 2 








K o ( i 


exp 


.2 V 







*2 



dsi / ds 2 



o 



T 



xF 



(gl ~ S2)Vq" 



.si- s 2 
2 +l % 



- s x )a{t 2 - s 2 ), 
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where a(t) = exp(iekt/h)Uk(t) (note a{t) is state independent), for which an equation 
similar to Eq. (5) can be derived. In this case, there is only one equation to be solved, 
instead of iV-coupled equations (cf. Eq. (7)), which results in a significant simplification 
of the entire formalism. Moreover, various analytic solutions can be obtained, as we 
exemplify below, by analyzing the adiabatic and the diabatic evolutions of the occupation 
numbers for a system with a realistic level density of the form p(e) = poexp(f3e). (3 = 
corresponds to the case we have just described, of constant average level density, while 
the case of finite (3 approximates fairly well a many-fermion system. 



The adiabatic limit corresponds to hqXq/%Vq S> 1 (and also 1), from which we 



find 



exp 



f 2irW 

{-IT 



rVo 
X n 



1 



t 



2ttW \t\ 

h 



(12) 



All odd moments of H vanish identically (since F(x) = F(—x)), and in the limit t — > oo, 
all even cumulants of Hq increase linearly in time. If F(x) = exp(— x 2 /2) (we shall use 
this form hereafter for illustrative purposes) then 

2n 



D(V 2 ) 



2nW n t fhVn 



(2n)\ 



x 



2 n n\ 



hw v 2 hr l v 2 



y 2 r, v2 ? 

(r^ = 2tcWq) resulting in a non-Gaussian distribution. 



(13) 
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In the diabatic limit, KoXo/hVo < 1, we find 



exp 



h 2 V 



IT 



exp 1 T I 2 7T 



exp 




In this case again all the cumulants of Ho increase linearly in time 

d2 2 \ / • \ n / • 



((mrnmt)}) 



h'Vo 6XP l 8 JL/5 



H, 



2y/2, 



t. 



(14) 



(15) 



where H n (x) are Hermite polynomials. From the explicit expressions for the first and 
second cumulants we thus obtain that 



(3D = 7 1 + 



(16) 



which shows that the Einstein fluctuation-dissipation theorem is generally violated. 

The most salient feature of the solutions (12) and (14) becomes evident when one 
considers the asymptotic behaviour of the cummulants. Since cumulants of higher than 
second order are nonvanishing, Gaussian processes are not obtained in any of these lim- 
iting cases for the energy diffusion. ((■?/>(£) ?/>(£))) either increase indefinitely with n 
or increase subsequently after an initial decrease, depending on the values of parameters. 
As a result the energy distribution has very long tails. In particular for Vq = the dis- 
tribution corresponding to Eq. (12) has a Lorentzian shape. These features imply that 
a Langevin or Fokker-Planck approach to energy dissipation is at least questionable. As 
we have discussed in Ref. QXUJ the present results apply equally to the GUE case. 

The present approach treats the fast subsystem quantum mechanically and the slow 
subsystem classically, as has been done often in the past [0 . The energy diffusion process 
is described in terms of some intrinsic characteristics of the many-body system (thermo- 
dynamic temperature, spreading width IN-, k and X ) and Vq. It is not clear yet whether 
these characteristics have a meaningful classical limit separately or only in a given com- 
bination, and this seemingly points to an apparent lack of a classical limit for the fast 
degrees of freedom (h — > 0) of the solutions (13), (15). In Ref. |J, friction was obtained 
only in a classical treatment of both fast and slow system, while in Ref. || , dissipation 
and friction appear only in a explicit quantum treatment (path integral) of the entire 
system and the presence of quantum fluctuations in the slow subsystem was essential. 
The wall formula || leads to a diffusion constant D oc V 2 as we have obtained here (see 
Fig. 3 and Eq. (13)), but is essentially a classical result, which does not depend in any 
significant way on h, and apparently reflects a different underlying mechanism. It will be 
highly desirable to identify the classical limit of the present approach. 

In conclusion, we have presented numerical and analytical solutions of the time de- 
pendent evolution of a driven complex quantum system, such as a nucleus, under the 
assumption that the number of levels is large. The parametric random matrix approach 
chosen here incorporates the essential attributes of the intrinsic dynamics, namely: an 
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exponentially increasing level density, GOE spectral fluctuations and loss of correlations 
during large amplitude collective motion. We have shown that the resulting energy dif- 
fusion process is highly non-Gaussian in character, that the energy distribution has very 
long tails and also that the energy diffusion constant is proportional to the square of the 
collective velocity. 
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Figure 1: Instantaneous eigenvalue spectrum E n (X) as a function of the "shape" (X), 
for a Hamiltonian of the ensemble defined by Eqs. (1-3), with [i?o]jfe — k5jk, using the 
correlator F(x) = exp(— x 2 /2) (top) and F(x) = exp(— \x\) (bottom). 
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Figure 2: The time dependence of the occupation probabilities nk{t), for k — 1, . . . , 51 
(in this case rik(t) = nio2-fc(0)> where k counts from bottom to top in the figure, for the 
case of fast, V = 4 (top), and slow, V = 1/16 (bottom), driving velocities. 
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Figure 3: The time dependence of the energy variance A E (t) for a range of velocities 
V = 4, 2, 1, 0.5, 0.25, 0.125 and 0.0625. The highest curve corresponds to the larger 
velocity. The insert shows the diffusion constant D(Vq) as a function of Vq, indicating 
D oc \/ 2 . 
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